The van der Waals energy of atomic systems near absorbing and dispersing bodies 
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Within the frame of macroscopic quantum electrodynamics in causal media, the van der Waals 
interaction between an atomic system and an arbitrary arrangement of dispersing and absorbing 
dielectric bodies including metals is studied. It is shown that the minimal-coupling scheme and the 
multipolar-coupling scheme lead to essentially the same formula for the van der Waals potential. 
As an application, the vdW potential of an atom in the presence of a sphere is derived. Closed 
expressions for the long-distance (retardation) and short-distance (non-retardation) limits are given, 
and the effect of material absorption is discussed. 
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I. INTRODUCTION 

The involvement of van der Waals (vdW) forces in a 
variety of physicochemical processes and promising po- 
tential applications (such as the construction of atomic- 
force microscopes pj or reflective atom-optical elements 
0) have created the need for a very detailed understand- 
ing and controlling them. Casimir and Polder [j| were the 
first to study the vdW interaction within the frame of 
rigorous quantum electrodynamics (QED). Investigating 
the interaction between an atom and a perfectly reflecting 
semi-infinite (planar) body, they found that in the short- 
distance (non-retardation) limit the interaction poten- 
tial U(z) behaves like z~ 3 (z, distance between the atom 
and the interface) , whereas in the long-distance (retarda- 
tion) limit it behaves like z~ A . In their theory, Casimir 
and Polder quantized the (transverse) vector potential 
of the electromagnetic field (outside the perfectly reflect- 
ing body) in terms of normal modes and coupled them 
to the atom according to the minimal-coupling Hamilto- 
nian, with its Coulomb part being determined by means 
of the method of image charges. They then calculated 
the (lowest-order) change of the ground-state energy of 
the system arising from this coupling, which is a func- 
tion of the atomic position and thus plays the role of the 
potential energy that determines the force acting on the 
atom. 

In the earlier experiments [4| , which studied the deflec- 
tion of thermal atomic beams by conducting surfaces, the 
observed signal was extremely low. Nevertheless, qualita- 
tive trends in agreement with the z~ 3 law were observed. 
Only recent progress in experimental techniques has ren- 
dered it possible to detect vdW forces with sufficiently 
high precision 

H IS II II In particular, by using 
atomic passage between two parallel plates, the z~ 4 re- 
tarded potential could be verified yj . Other methods for 
measuring vdW forces have been based on transmission 
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grating diffraction of molecular beams |6g , atomic quan- 
tum reflection Q, evanescent-wave atomic mirror tech- 
niques 181, and (indirect) measurements via spectroscopic 
means [9|. Proposals have been made on improvements 
of monitorin g th e vdW interaction by using atomic in- 
terferometry (1C|. 

Since the appearance of Casimir's and Polder's pio- 
neering article in 1948 there has been a lar ge b od y of 
work on the vdW interaction (see, e.g., jllL Ha . Ha . HH llq] 
and references therein). Roughly speaking, there have 
been two routes to treat the problem. In the first, 
which closely follows the ideas of Casimir and Polder, 
explicit field quantization is performed by applying stan- 
dard concepts of QED, such as normal-mode techniques 
[S [H E d M I2l|j2pj23l| . In particular, extensions 
of the theory to one [171 Hq| and two semi- infinite dielec- 
tric walls [l3.l20| . thin metallic films |2l|, and cylindrical 
and spherical dielectric bodies [2^ have been given, and 
the problem of force fluctuations on short time scales 
has been studied j2^. The calculations have typically 
been based on macroscopic QED, by applying normal- 
mode decomposition and including in it the bodies by 
the well-known conditions of continuity at the surfaces 
of discontinuity. Since in such an approach the depen- 
dence on frequency of the response to the field of the 
bodies cannot be properly taken into account, material 
dispersion and absorption are ignored. The problem does 
of course not occur in microscopic QED, where the bod- 
ies are treated on a microscopical level by adopting, e.g., 
harmonic-oscillator models (see, e.g., |l7|). Apart from 
the fact that the calculations are rather involved, the re- 
sults are model-dependent. 

To overcome the difficulties mentioned, in the second 
route, the calculations are based on linear response the- 
ory, without (exp licitly) quantizing the electromagnetic 
field |H )M El Hi |H H 111. All the relevant 
entities are expressed in terms of correlation functions 
which in turn are related, via the fluctuation-dissipation 
theorem, to response functions. The method has been 
employed to investigate the vdW energy of an atom near 
a semi-infinite (planar) body made of dielectric E3 ; 
metallic [2M [23, ionic-crystallic [2^, and birefringent 
material |29j. Further, the problem of an atom near a a 
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dielectric plate 26] , a metallic sphere 30] , and a nanowire 
pHf has been considered. Since the method borrows from 
equilibrium statistical mechanics, its applicability is re- 
stricted to the vacuum and thermal quantum states. 

In this article, we give a unified QED approach to 
the vdW interaction between an atomic system (such 
as an atom or a molecule) and dispersing and absorb- 
ing dielectric bodies including metals. Starting from the 
quantized version of the macroscopic Maxwell field, with 
the medium being described in terms of a spatially vary- 
ing, Kramers-Kronig consistent (complex) permittivity, 
we derive an expression for the vdW potential that ap- 
plies to arbitrary body configurations. The formalism 
can be regarded as being a generalization of the normal- 
mode formalism of macroscopic QED, so that it can be 
applied to other than the vacuum and thermal states, and 
it also allows extensions like the inclusion of dynamical 
interactions between the atomic system and the medium- 
assisted electromagnetic field. Roughly speaking, the 
mode expansion is replaced with some source-quantity 
representation in terms of the Green tensor of the macro- 
scopic Maxwell equations, in which material dispersion 
and absorption are automatically included. Further, we 
close the gap between the minimal-coupling scheme and 
the multipolar-coupling scheme by showing that both ap- 
proaches lead to equivalent results. 

To give an application, we consider the vdW poten- 
tial of an atom in the vicinity of a dispersing and ab- 
sorbing microsphere. Microspheres may be interesting 
candidates for QED experiments (see, e.g., [32l l33l] and 
references therein) with atomic beams. The atom-surface 
distances should be adjusted so that the atoms fly close 
enough to the surface to facilitate a strong coupling with 
the microsphere resonances, but not get adsorbed on the 
microsphere surface because of the attractive vdW force. 

The article is organized as follows. In Section [H] the 
formalism is outlined and an expression for the vdW po- 
tential is derived. In Section IlIII the formalism is applied 
to an atom near a sphere. A summary and some conclu- 
sions are given in Section llVI 



II. THE VAN DER WAALS ENERGY 



A. The quantization scheme 



Let us consider an atomic system (such as an atom or a 
molecule) interacting with the quantized electromagnetic 
field in the presence of macroscopic, dispersing and ab- 
sorbing dielectric bodies. In the nonrelativistic limit, the 
minimal-coupling Hamiltonian in Coulomb gauge reads 



34, 3f 
H = 



dujhu J d 3 r&(r,w)f(r,u) 



E 



2m r 



p a - q a A(r A + r a ) 



+ | / d 3 rp A (r)<p A (r) + J d 3 r^(r)^ M (r), (1) 

where ra a and q a are respectively the masses and charges 
of the particles constituting the atomic system, while f Q 
and p a are respectively their coordinates (relative to the 
center of mass r A ) and canonically conjugated momenta. 
The first term in the Hamiltonian describes the combined 
system of the electromagnetic field plus the macroscopic 
bodies (including dissipative systems) in terms of bosonic 
vector fields f(r, u>), which satisfy the commutation rela- 
tions 

7i(r,w),/f(r , J w , )l =<M(r -r')5(u;- a/), (2) 



/ i (r,«),/ i (r',w / 



/7(r, W ),£(rV) 



= 0. 



(3) 



The second term is the kinetic energy of the charged 
particles, while the third term describes their mutual 
Coulomb interaction, where 



and 



<^A(r) 



i 



47T£ 



,3 i PA(r') 
d r -T 



q a S[r - (r A + r a )] 



(4) 



(5) 



are respectively the scalar potential and the charge den- 
sity of the atomic system. The last term accounts for the 
Coulomb interaction of the particles with the medium. 
The vector potential A(r) and the scalar potential (/3m ( r ) 
of the medium-assisted electromagnetic field are given by 



A(r) = / dw(io;) _1 E- L (r,a;)+H.c. ) 



V<p M (r)= 6ljW(t,u}) +H.c. 



(6) 



(7) 



where the symbols _L and || are used to distinguish trans- 
verse and longitudinal vector fields, respectively. In par- 
ticular, E^(r, lj) and E"(r, oS) read 



& ± W(r,w)=J dV <5 ±(ll) (r-r')E(r», 



(8) 



where <5 X (r) and <5"(r) are the transverse and longitu- 
dinal dyadic 8- functions, respectively, and the medium- 
assisted electric field in the w-domain, E(r, to), is ex- 
pressed in terms of the basic variables i(r,u>) as 

E(r,w) = iiL*-" f dV \Ame(r',w) 

V tt£o c 2 J ( 9 ) 

x G(r,r»f(r» 
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(c= l/y/eono). Here, G(r, r',ui) is the classical Green 
tensor and Ime(r, uj) is the imaginary part of the com- 
plex, space- and frequency-dependent (relative) permit- 
tivity 



e(r, uj) — Ree(r, uj) + ilme(r, uj). 



(10) 



The Green tensor, which obeys the inhomogeneous par- 
tial differential equation 



For the following it is convenient to decompose the 
Hamiltonian Q as 



where 



H = H F + H A + H AF , 



(21) 



H F = duhioj d 3 r&(r,uj)f(r,u;), (22) 



V x V x — - e(r, uj) 



G(r,r',uj) = 8(r-r') (11) 



together with the boundary condition at infinity, has the 
following general properties (see, e.g., [35|): 



G*(r,r» = G(r,r',-^), 



(12) 



^E^ + Hd 3 r/3A(r)^ A (r) 



2m a 
f. 2 

y v^ + y 



2m a ' Aneo \r a - rp\ ' 



(23) 



Gij(r, r',u) = Gji(r',T,u), 



(13) 



^ J d 3 sei(s,w)G(r,s,w)G*(s,r',w) = ImG(r,r». 

(14) 

In this way, all electromagnetic-field quantities can be 
expressed in terms of the fundamental fields f(r, u>). In 
particular, the operator for the electric field reads 

E(r) = -l[A(r),£] -V0M(r)-V0 A (r) 

= E M r+> 1 ■ - —r , (15) 

^ 47re r- (r A + r a )l d 



where 



E M (r) = / dwE(r,u)+H.c. 
Jo 



(16) 



and the polarization field associated with the dielectric 
medium is given by 



where 



P M (r) = / dwP(r,w)+H.c, (17) 
Jo 



P(r, uj) = X (r, w)e E(r, uj) + P N (r, uj). (18) 



H 



AF 



d 3 r p A (r)0 M (r) - V — A(r A + r a ) 
p Q - ig«A(r A +f a )l . (24) 



Obviously, 77p is the Hamiltonian of the medium-assisted 
electromagnetic field, Ha is the Hamiltonian of the 
atomic system with eigenstates |n) and eigenvalues E n 
according to 



H A \n) = E n \n), 



(25) 



and Haf is the interaction energy between them. For a 
neutral atomic system in the electric-dipole approxima- 
tion, the latter simplifies to 



tl AY — tl K¥ + H AF , 



(26) 



Hll = -^^p Q A(r A )+dV^ M (r A ), (27) 



(28) 



where d is the electric dipole operator of the atomic sys- 
tem 



d = ^2q a r a . 



(29) 



Here, 

X(r,uj) = e(r,uj) - 1 
is the dielectric susceptibility and 



(19) 



P N (r,w) = *y^Ime(r,w)f(r,a;) (20) 
is the so-called noise polarization. 



B. The vdW energy in the minimal-coupling 
scheme 



Following the original line of Casimir and Polder y], 
we calculate the vdW energy of an atomic system as the 
position-dependent part of the leading-order correction 
to the unperturbed ground state energy due to the per- 
turbation according to Eq. I|26|l [together with Eqs. I|27|) 
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and (|28|l ]. Since the diagonal matrix elements of H^p, 
Eq. (|27|) . are zero, the non- vanishing contribution to the 
energy correction from first-order perturbation theory is 
due to Hj$, Eq. 1(25)1, Inspection of Eq. (J2HI shows 
that this contribution is of linear order in the fine struc- 
ture constant a = e 2 /(2eohc). We have thus to include 
that contribution from second-order perturbation theory 
which is of the same order in a. Therefore, we apply 
first-order perturbation theory for H^p and second-order 

perturbation theory for -H^f- The energy correction to 
the ground state thus reads 



where 



AE ~ + A 2 E, 



A,E ee (0|({0}|ff^|{0})|0) 

= (0|({0}|^^A 2 (r A )|{0})|0) (31) 



(30) 



and 



l3f |(0|({0}|if«|{l(r,a;)})|n) 



A 2 E = V / du f d 3 . 

h ^ i~ LU n + LU J ^ m a 



w n 



■dV^ M (r A )|{l(r,w)})|n) 



(32) 



Here, |0) and |{0}) are respectively the ground state of 
the atomic system and the ground state of the medium- 
assisted electromagnetic field, and 



\{1(t,u)})=F(t,w)\{0}) 



(33) 



denotes single-quantum Fock states of the fundamental 
fields. Further, 



w„ = (E n - E )/H 



(34) 



are the transition frequencies between the excited atomic 
states and the ground state. Note that due to the linear 
dependence of the vector potential and the gradient of 
the scalar potential on f(r,cj) [and V(r, u>)], the nonva- 
nishing matrix elements of H^p in the Fock-state basis 

(which is defined as the set of eigenvectors of Hp) are 
those between Fock states that just differ in one quan- 
tum. 

Expressing in Eq. I|31|l A(r A ) in terms of f(r) [and 
f^r)], on using Eqs. ©, ©, and ©, recalling Eq. (J33), 
and exploiting the integral relation l|14|) for the contrac- 
tion of two Green tensors, we derive after some lengthy 
but straightforward calculation 

AiE =^J2JL f°d W Im^(r A ,r A , W ), (35) 
7T ^ 2m a J 



where 

J-(II)G-L(ll)( rj r» 

= /d 3 s| dV ^ (ll) (r-s)G(s,s',^)<5 ± (H)(s'-r'). 

(36) 

In Eq. I|35J) and below, summation over repeated vector 
indices is understood. Using the sum rule 



'5] 2m~^ = 2fi E/ U n{don,id n O,j + do n ,jd n o t i) (37) 



2h 



(for a proof, see Appendix lA"|l . where 

do n ,i = (0\di\n), 



(38) 



denote the matrix elements of the electric-dipole operator 
(|29|l . we may equivalently represent Eq. (|35|> in the form 
of 

Al £=^°y;A % rdcIm^(r A ,r A , W ) 

„p roc 

= — VW dww„d 0n Im- L G ± (r A ,r A ,w)d„o- (39) 
In a similar fashion, we find that Eq. (|32|) leads to 



AoE ■■ 



Mo 

7T 



E 







/ da; 
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_0J n + W J 



d 3 s 



x y d 3 s' d „/x(s, w)Im G(s, s', w)/i(s', w)d„ 



(40) 



where, on using the relation l|A3(l . the matrix elements 
of the electric-dipole moment l(5S)l have been introduced 
and the abbreviating notation 



/j,(t,w) = S ± (r - r A ) 



■6H 



r - r A l 



(41) 



has been used. 

We substitute Eqs. (O and gDJ into Eq. l(3UJl. and 
obtain, on recalling Eq. (|41(l . 



AS 



7T 



y-v /" — ^_do„{w 2 Im ll G ll (r A ,r A ,o;) 
„ Jo u n + u L 

-w„w[lm- L G- L (r A ,r A ,Lj) + Im^G 11 (r A , r A , w) 
+ Im ll G ± (r A ,r A ,o;)]}d n o. (42) 



Let 1Z be the (small) region of space where the atom is 
situated and the permittivity can be regarded as being 
effectively not varying with space, i.e., e(r, ui) ~ e(r A , ui) 
if r £ 7Z. For r, r' S 1Z, the Green tensor can then be given 
in the form 

G(r, r', u) = G^ (r, r', w) + G«(r, r', «), (43) 



5 



with (r, r', to) denoting the (translationally invari- 
ant) bulk Green tensor that corresponds to e(T^ : to) and 
G^'(t,t',lo) being the scattering Green tensor that ac- 
counts for the spatial variation of the permittivity. In 
practice, the atom is typically situated in a free-space 
region, so that G^ is simply the vacuum Green tensor. 
According to the decomposition of the Green tensor in 
Eq. ((43(1 . the energy correction AE given by Eq. 142|) 
consists of two terms, 

AE = A£(°) + AE^(r A ), (44) 

where the rA-fndependent term AE^°\ which is related 
to the bulk Green tensor, gives rise to the (vacuum) Lamb 
shift, whereas the rA-dependent term AE 1 - 1 ") (r A ), which 
is related to the scattering Green tensor, is just the vdW 
energy sought: 

U(v A ) e, A£«(r A ) =-»Vf d „ 

n „ Jo w„ + w 



To further evaluate this expression, it is convenient to 
express it in terms of the whole scattering Green ten- 
sor rather than its imaginary part. For this purpose, 
we write ImG (1) = (G$ - G (1) *)/(2i), recall the rela- 
tion (|12|l . and change the integration variable from — lo 
to lo. Equation (|45|) then changes to 



x |w 2 WgW (rA, r A , lo) - lo u oj 

x [im^G^^^A, r A , to) + Im X G (1 >H (r A , r A , w) 

+ ImllG« ± (r A ,r A ,cj)]}d n o. (45) 



U(ta) 



P- Vd J r -^-\^G^\v A ,v A ^) + ^G™{r A ,v A ,u) + WG^{v a ,v a ,lo)} 

2lTT ^ \J 10 n + LO I 

-Lo^G^Hr A ,r A7 co)}+ f -^L n ^ G ^(r A ,r A ,u) + ^G^(r A ,r A , to) + UG«> A ,r A , lo)] 

> J -co U n - LO I 

+w 2|l G (1)ll (r A ,r A) w)}]d n o. (46) 



This equation can be greatly simplified by using contour- 
integral techniques. Note that G l -- 1 \r A ,r A ,Lo) is an ana- 
lytic function in the upper complex half plane (Imw > 0, 
lo^O). Further, ± GW x (r A , r A , lo), x G^W(v a ,t a ,lo), 
and llGrW- L (r A , r A , lo) tend to zero as lo approaches zero, 
because the asymptote of the Green tensor contains no 
transverse components [cf. Eq. (|B4|) ]. Finally, the term 
lo 2 HG' 1 ^l(r A ,r A ,w) is also well-behaved for vanishing lo, 
as can be seen from Eq. (|B7ll . Consequently, the inte- 
grands of the w-integrals in Eq. 1)46(1 are analytic func- 
tions without poles in the whole upper complex half 
plane, including the real axis. We may therefore apply 
Cauchy's theorem, and replace the integral over the pos- 
itive (negative) real half axis by a contour integral along 
the positive imaginary half axis (introducing the purely 
imaginary coordinate lo = iu) and along a quarter circle 
with infinite radius in the first (second) quadrant of the 
complex frequency plane. Since the integrals along the 



infinitely large quarter circles vanish [cf. Eq. I|B3|) ]. we 
finally arrive at 



U(r A ) 



Mo 



E 



du 



rd „G (1) (r A , r A , iu)d n0 , 
(47) 

where the identity G^ = ± G^- L + ^G™ + \\G^ 1]± + 
"G^ 1 )" has been taken into account. 

Introducing the (lowest-order) ground-state polariz- 
ability tensor 



a(uo) = lim — ■ 



irjio 



■ d „ ® d„ (48) 



of the atomic system (see, e.g., 36]), we may represent 
Eq. 147(1 in the equivalent form of 

U(r A ) = ^ r duu 2 a l0 {iu)Gfhv A ,v A ,iu). (49) 
27T Jo 
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It is worth noting that Eq. i|49|) directly follows from QED 
in causal media, without the need of additional assump- 
tions borrowed from other fields. Equation 1(49(1 expresses 
the vdW potential of an arbitrary atomic system (such 
as an atom or a molecule) in the presence of an arbitrary 
configuration of dispersing and absorbing macroscopic di- 
electric bodies in terms of the polarizability tensor of the 
atomic system in lowest order of perturbation theory and 
the scattering Green tensor of the macroscopic Maxwell 
equations. 

In particular for an atom, one can make use of the 
spherical symmetry and reduce Eq. ((49(1 to 



U(ta) = ^7 J duu 2 a(m)G l ( l 1) (rA,r A ,m), 



(50) 



where 



a(uj) 



lim — 



£ 



uj* — irjui 



l 0n\ 



(51) 



This result agrees with the results inferred from (se- 
mi-) classical linear response theory. Note that the field 
susceptibility introduced in Refs. |2J] and |2fj differs from 
the scattering Green tensor by a factor of J 1 . Needless to 
say, that in the special case of a two-level atom, Eq. 1(50(1 
reduces to the result, e.g., in Ref. |25j. The derivation 
of Eq. I|5U|) shows that it can be regarded as the natu- 
ral extension of the QED results obtained on the basis 
of the normal-mode formalism, which ignores material 
absorption. 



being the polarization associated with the (neutral) 
atomic system. Using H from Eq. JQ), we derive [35ll37| 

H= [d 3 r [ duo huP (r,Lj)f(r,ij) + Y 7T—{p a 
J Jo „ 2m a I 



dAAr a x B [r A + Af Q ] 



d 3 rP A (r)P A (r)- d 3 rP A (r)E„(r), (55) 



25, 



where B(r) = Vx A(r) with A(r) from Eq. © [together 
with Eqs. © and ®], and E M (r) is defined by Eq. JTSJ) 
[together with Eq. JJ]J]. Note that in the multipolar- 
coupling scheme the operator of the electric field strength 
is defined according to 



Efr) 



1 

ih 



A(r),H - V£ M (r) - V<£ A (r) 



1 



E M (r) P A (r), 



£ E M (r) =e E(r) + P A (r). 



(56) 



(57) 



Hence, £oE]v[( r ) has the meaning of the displacement field 
with respect to the polarization of the atomic system. 

In the electric-dipole approximation, Eq. 1)55(1 simpli- 
fies to 



Tt = Hf + Ti-A + ~Hkf , 



(58) 



The vdW energy in the multipolar-coupling 
scheme 



Let us turn to the multipolar-coupling scheme widely 
used for studying the interaction of electromagnetic fields 
with atoms and molecules. Just as in standard QED, so 
in the present formalism the multipolar-coupling Hamil- 
tonian can be obtained from the minimal-coupling Hamil- 
tonian by means of a Power-Zienau transformation, 



H = U^HU, 



where 



U 



exp 



d 3 rP A (r)A(r) 



(52) 



(53) 



with 



q a r a / dAJ[r-(r A + Af 



(54) 



where 



and 



H F = Jd 3 rJ^ dw/ia>ft(r,cj)f(r,u;) 
*A = £|t + ^ /^P A (r)P A (r), 



(59) 



(60) 



respectively, are the unperturbed Hamiltonians of the 
medium-assisted electromagnetic field and the atomic 
system, and 



Haf = -dE M (r A ) 



(61) 



is the interaction energy between them, where d is the 
atomic dipole operator given by Eq. Recall that 

Em(i" A ) must be thought of a being expressed in terms of 

the fundamental field variables f(r a) [andF(r A )]. Com- 
paring Eq. l(6Tj|) with Eq. iffify . we see that, on taking into 
account the the relationship 



4 /d 3 ^ A (r)^ A (r) = l / d 3 rp[(riP^r i. ((,2) 
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the atomic Hamiltonians Tl A and H A are different from 
each other, so that the solution of the eigenvalue problem 



H A \n') = E'Jn') 



(63) 



may be different from that defined by Eq. (|25|l . Keeping 
in mind this difference, we drop the primes denoting the 
atomic eigenvalues and eigenstates from here on. 

In contrast to the interaction energy in the minimal- 
coupling scheme, Eq. (|26() . the interaction energy in the 
multipolar-coupling scheme, Eq. (|61|) . is linear in f(r A ) 
and f f (r A ). As a consequence of the latter, there is no 
first-order correction to the ground state energy. We thus 
have 



where 



A?E 



E 



AE ~ A 7 E, 



dtu / d r 



(64) 



|(0|({0}|dE M (r)|{l(r,c)})|n) 

E - (E n + hw) 



(65) 



[cf. the first line in Eq. I1.S2I) with Haf instead of H^L]. 
In complete analogy to the derivation of Eq. (|4t)|l , we find 
that 



A 2 E 



Mo 

7T 



E 



dtu ■ 



— do„ImG(r A ,r A ,w)d„o ■ 

- U! 

(66) 

To further evaluate the rA-dependent part U(r A ) = 
A^ E(r A ) of A?E, which results from the scattering part 
of the Green tensor, G^^a, r A , w), and gives the vdW 
energy, we again write ImG' 1 ' = (G*- 1 -' — G^ 1 '*)/(2i), use 
the relation (|12|) . and change the integration variable 
from —uj to lo. After some algebra we arrive at 



U(r A ) 



Mo 

2l7T 



E d °« / 



duj- 



G«( 



r A ,r A ,uj) 



o 



duj- 



G«( 



r A ,r A ,o;j 



-InO 



(67) 



As we already know, the integrands of the two frequency 
integrals appearing in Eq. 1671) are analytic functions in 
the upper half of the complex frequency plane, including 
the real axis [cf. Eq. JHJJ]. We therefore can apply 
contour integral techniques in a similar way as in the 
derivation of Eq. lj4*7|) from Eq. lj4U|) . It is not difficult to 
see that the result reads 

f7(r A ) = — V* / du UJnU d 0ra G (1) (rA,rA,m)d w0 , 
^ „ Jo ul+u 2 

which has exactly the same form as the minimal-coupl- 
ing result (|47|) . so that it can also be given in the form 



of Eq. H49|) . Recall that the values of u n and do„ ob- 
tained in the minimal-coupling scheme may be different 
from those obtained in the multipolar-coupling scheme, 
because of the somewhat different eigenvalue equations 
(gSJ) and (EH). 



III. APPLICATION: AN ATOM NEAR A 
SPHERE 

Let us apply the theory to an atom near a dispers- 
ing and absorbing dielectric (micro-)sphere surrounded 
by vacuum. The material of the sphere of radius R is 
assumed to be homogeneous and isotropic, having a per- 
mittivity s(oj) . The coordinate system is chosen such that 
its origin lies at the center of the sphere. The scattering 
Green tensor can be given by [3g| 



(n — to)! 
(n + to)! 



Attc ^ ^ K n(n + 1) 

n=lm=0 v ' 



Bf} ^2 M„ m ,p(r A ) ® M 

p=— 1,1 



nm.p I TA J 

p=-l,l 



(69) 



where M„ mjP (rA) and N nmtP (r A ) are even (p—1) and 
odd (p = — 1) spherical wave vector functions, which can 
be expressed in terms of spherical Hankel functions of 
the first kind, hn\r), and associated Legendre functions, 
P™(cos9), in the following way: 

M nm , ±1 (r A ) = T-^tjt h£\k r A )P™(cose) 

sin , .as,, s dP™(cos6») cos, , s 
x cos (™4>)ee ~ h^{k r A ) ^— (to^, 



N nm , ±1 (r A ) = n(n + 1) h ^} k ° r ^ P ™ icost 



(70) 



x . lm<p)e r 
sm v ^' 



kor A 

1 dfrA^CAorA)] 



k r A 



dr A 



dP™(cos( 



x . (m(j))eg =f P" 1 (cos 9) (mdfiej, 
sin v ^' sin# cos v ' v 



(71) 



Here, ko=Lu/c = iu/c is the vacuum wave number, and 
e. r , eg, e^j are the mutually orthogonal unit vectors 
pointing in the directions of r, 9, and </>, respectively. 
The coefficients B^f and in Eq. ^ read 



[Zljn(zi)] 'jn(zp) ~ [zojnjzo)]' jnjzi) 

[zijn(zi)]' h£ \z ) - [z h£ \z )]'j n (zi) 



, (72) 



8 



B 



N 



B»{iu) 

e(iu)j n (zi) [z j n (zp )] ' - j n (z ) [zij n (zi)]' 
e(iu)j„(zi)[zoh£\zo)]' ~ h£'(z ) [zij n (zi)]' 



(73) 



where z = k R and z\ = kR, with k — k ^e(iu) being the 
wave number inside the sphere, and j n (z) is the spherical 
Bessel function of the first kind. The primes indicate 
differentiations with respect to Zq or Z\ , respectively. The 
coefficients represent contributions from transverse 
electric (TE) waves reflected at the surface of the sphere, 
while the coefficients B^ represent those from transverse 
magnetic (TM) waves. 

Substituting the trace of G^\rx, r A , to) from Eq. i|tj9|) 
[together with Eqs. (|7UJ) and f7T)l] into Eq. (JHOJ) yields the 
vdW energy sought. The sums over p can then easily be 
performed using the orthogonality of the unit vectors e r , 
eg, and e^, and the sum over m can be performed with 
the aid of the summation formulas in Appendix [U] So 
after a lengthy, but straightforward calculation we arrive 
at the following result: 



U(r A ) = 



r du u 3 a(iu)Y(2n + l) 



x <^ B 



M 



Of A; 



+ n{n + 1)B 



N 



hn\k r A ) 



kor A 



+ B 



N 



1 d[r A h { n\k r A )} 



k r A 



dr A 



(74) 



In this case, Eq. I|74|l reduces to 
hcR 3 1 



U(r A ) 



4ir 2 e rl 



A -'o 



dz a{icz/r A ) 



e(icz/r A ) — 1 
s(icz/r A ) + 2 



2(1 + zf + (1 + z + z 2 f 



(76) 



(see Appendix lD|l . where it turns out that the TE waves 
do not contribute. Since the inequalities e(icz/r A )>l 
and a{icz/rx) > are valid, the vdW potential is nega- 
tive, and the resulting force between the atom and the 
sphere is attractive. 

As is seen from Eq. ()76|l , the main contribution to the 
integral comes from the region where z < 1/2. There- 
fore, for sufficiently large distances, the contributions 
from small frequencies dominate, and we can (approx- 
imately) replace the atomic polarizability and the ma- 
terial permittivity in Eq. (|76|l with their static values 



a(cu = 0) and e^ ) = e(uj — 0), respectively. The in- 



tegration can then be performed in closed form to yield 
the asymptotic distance law 



U(r A ) = - 



167T 2 £ 



e(°) + 2 rl 



in this so-called retarded limit. Note that in the opposite 
nonretarded limit, where the contributions of a(to) and 
e(ui) at all frequencies have to be retained, Eq. (|7B^I re- 
duces to the result given in Ref. j^, where ar^ 6 law was 
found. The (formal) limit e^ ) -^oo in Eq. fTTjl obviously 
corresponds to a metallic sphere 



U(r A ) 



23hcR 3 aW 1 
167r 2 £o rl 



(T_A_ 

\ R 



(78) 



Note that the vdW potential does not depend on the an- 
gle variables of the atomic position, but only on the dis- 
tance of the atom from the center of the sphere, as can 
be anticipated from the symmetry of the system. Re- 



call that the terms proportional to B, 



M(N) 



represent the 



contributions from the TE (TM) waves. Equation 174|) 
applies to an arbitrary dielectric sphere. In particular, 
when material absorption is omitted, then the result in 
Ref. [Hi can be recovered. 



1. Long-distance limit 

A detailed analysis of Eq. l|74Jl requires numerical com- 
putation. Here, however, we would like to focus our at- 
tention on two interesting limiting cases, where the atom 
is very far from or very close to the sphere. Let us first 
consider the limit of the atom being far away from the 
sphere, 



r A > R. 



(75) 



Note that the decrease of the force with the distance is of 
three powers stronger than in the case of the atom being 
near a planar body. 

In particular, if we introduce the static polarizability 
of the sphere (see, e.g., [3^1 



(0) A £ 



(0) 



1 



(0) 



R 3 



we may rewrite Eq. 177(1 as 



U(r A ) 



a (0 U (0) 

sph 
(47T£ ) 



23hc 1 



2 47r r\ 



r_A_ 
R 



(79) 



(80) 



Interestingly, Eq. I|80|) also applies to the vdW potential 
between two atoms 3] , if the (static) polarizability of the 
sphere is replaced with the polarizability of the second 
atom. 



2. Short- distance limit 

Let us now proceed to the short-distance limit of the 
atom being located at a position very close to the sphere, 
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R 



< 1 



(81) 



(ArA = ta — R)- In this case, from Eq. l(7^|) it follows 
that 



U{ta) 



1 



167r 2 £ (Ar A ) 3 J 



' dua(iu) 6( ? U l 1 (82) 
v ; e(iu) + 1 



(see Appendix [EJl . Note that again the TE waves do not 
contribute to [/(ta)- 

As expected, the dependence on distance of the (at- 
tractive) vdW potential corresponds to that obtained in 
the case of the atom being near a planar body. In fact, 
it exactly looks like that derived in Ref. [2(| for an atom 
in the vicinity of a planar, semi- infinite, non-absorbing 
dielectric. In particular, if the (model) assumption [e(iu) 
— l]/[e(m) + 1] = 1 were made for all values of u, then 
Eq. JB2J would lead, on using Eq. (j5l")) . to the result Q 



U(r A ) 



(0|d 2 |0) 1 
48^£o (Ar A ) 3 



(83) 



With regard to Eq. I|82|) . we therefore find that, on re- 
calling that a(iu) > 0, 



d_ 



dU(r A ) 



dr A 



< 0. 



Hence, the vdW force monotonically decreases with in- 
creasing absorption constants. 




0.01 



0.015 



0.02 



3. Material absorption 

To explore the effect of material absorption, we may 
assume a permittivity of Drude-Lorentz type, 



e(w) = 1 + 



(84) 



FIG. 1: Absolute value of the normalized van der Waals force 
C\dU(r A )/dr A \ [C* = 167r 2 eo/(EJdon| 2 Af), \ t =2ivc/Qi] as 
a function of the atom-surface distance for various strengths 
of material absorption. In the calculation, a metallic per- 
mittivity of Drude type according to Eq. 1841 and a (de- 
generate) single-resonance atomic polarizability are assumed 
with oji = 0, u n /Q.i — 7x 10" 1 , and 7i/f2; = 10~ 2 (solid line), 
7z/fi; = 10 _1 (dashed line), and ^i/fli = 1 (dotted line). 



where u>i and 7/ are respectively the (transverse) res- 
onance frequencies and the associated absorption con- 
stants, and the frequencies tti are proportional to the 
so-called oscillator strengths. From Eq. (|84|l it is seen 
that in the limit lu — > the resulting static permittivity 



,(0) 



E 



2 



1 "1 



(85) 



is independent of the absorption parameters. Since it is 
the static permittivity that enters Eq. (|77() . we see that 
the long-distance asymptote of the vdW potential is not 
influenced by material absorption. 

With decreasing distance the range of frequency that 
must be taken into account increases. Thus, the fre- 
quency response of the permittivity becomes crucial to 
the strength of the vdW force. Let us consider the short- 
distance law 1 (52) 1 . From Eq. ((51)l it follows that 



d e(iu) — 1 
dji e(iu) + 1 

that is to say, 



1 



[e(iu) + l] 2 (uf + u 2 + uji)" 2 



d e{iu) — 1 
dji e(iu) + 1 



< if u > 0. 



(86) 



(87) 



Figurenillustrates the influence of material absorption 
on the vdW force acting on an atom located near a metal- 
lic sphere in the short-distance limit. It can be seen that 
the effect of material absorption increases with decreas- 
ing atom-surface distance. In particular, at a distance 
of ArA — 10 _2 A/, an increase of the relative absorption 
parameter from ji/fli = 10 -2 to J1/Q1 = 1 would reduce 
the magnitude of the force by nearly thirty percents. 



IV. CONCLUSIONS 

Within the frame of macroscopic QED, we have de- 
rived an expression for the vdW potential of an atomic 
system near an arbitrary configuration of dispersing and 
absorbing bodies. It generalizes the results obtained by 
means of normal-mode expansion and may be regarded 
as a foundation of the results inferred from linear re- 
sponse theory. We have performed the calculations for 
both the minimal-coupling scheme and the multipolar- 
coupling scheme and shown that the results are essen- 
tially the same. 

We have applied the theory to the vdW interaction 
between an atom and a sphere. From the integral ex- 
pression, we have derived the correct long-distance law 
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corresponding to the retardation limit and recovered the 
short-distance law corresponding to the non-retardation 
limit. In particular, replacing in the long-distance law 
the polarizability of the sphere with that of an atom just 
yields the vdW potential between two atoms. On the 
other hand, for sufficiently small distances of the atom 
from the sphere the vdW potential approaches the po- 
tential observed for an atom near a planar body. 

It is worth noting that in the long-distance limit it 
is the static permittivity that enters the vdW potential. 
Hence material absorption has no effect on it. However, 
with decreasing distance of the atom from the sphere the 
relevant frequencies extend for a finite (increasing) inter- 
val and material absorption becomes substantial, thereby 
diminishing the strength of the force. 

In this article, we have restricted our attention to 
ground-state systems and calculated the vdW potential 
in lowest-order of perturbation theory with respect to 
the interaction of the atomic system with the medium- 
assisted electromagnetic field. The theory allows of 
course extensions in several respects. As a consequence of 
the lowest-order perturbation theory, the energy denom- 
inators that enter the polarizability of the atomic system 
are the unperturbed ones, without consideration of the 
level shift and broadening caused by the presence of the 
bodies. In fact, the polarizability of an atomic system 
is expected to drastically change when it becomes close 
to a macroscopic body and the spontaneous decay thus 
becomes purely radiationlcss, with the decay rate being 
proportional to Ar^ 3 J22|. Since the level broadening 
is essentially determined by the spontaneous-decay rate, 
the polarizability becomes distance-dependent - an effect 
that needs careful consideration. 

Since the electromagnetic field in (linear) magnetic me- 
dia can be quantized analogously , another interesting 
extension of the theory be the inclusion in it of composite 
materials characterized by both a complex permittivity 
and a complex permeability. Interestingly, such materi- 
als, which have been fabricated recently, are left-handed. 
Last not least the underlying quantization scheme ren- 
ders it also possible to extend the theory to atoms and 
molecules in excited states and treat the motion of driven 
atomic systems. 
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APPENDIX A: DERIVATION OF EQ. STfh 



From the commutation relation 



Ma 



2m/3 



p a>i (Al) 



together with the eigenvalue equation 125fl we find that 

(0\p a ,i\n) = - l -^(0\[r a , h H A ]\n) 

= ^{Q\r ai H A - H A r a ,i\n) 

n 

= T^{E n -E o )(0\r a ^\n) 



(A2) 



Thus 



2h x — ' m a 

a n 



2h ^ m n 



E — (0|Pa,i|n) = -iu n q a {0\f aA \n) 

a a 

= -iu n (0\di\n). (A3) 
Using Eq. (|A3|I . we derive 

^E>» ((0|d>)(n|d>) + (0|d,|n)<n|*|0) 

n 

-<0|d>)(r#a,i|0) 

[Pa,i,dj]\Q) 
a 

-Eslb*" (A4) 

a 

which is just Eq. (|37|l . 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF 
THE GREEN TENSOR 

The asymptotic behavior of the Green tensor for large 
frequencies reads [3j| 

lim ^G(r,r»=-5(r-r'), (Bl) 



lim ^G (0) (r,r',w) = -tf(r-r'). (B2) 

\u)\ — >oo C 

If r and r' lie in a common region of constant permittiv- 
ity, we can use Eq. I)43|l . and subtract the two equations 
(|BT)) and (|B2|) to obtain 



lim ^-G (1) (r,r',w) = 0. 



\u!\—fOO C 



(B3) 
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In the low-frequency limit we have [3£ 

,2 



lim — G(r,r',w) = - H i" 1 H (r, r') , (B4) 



M-»o c 



where 



L(r,r') = lim / d 3 s <5" (r - s)e(s, w)<S" (s - r'). (B5) 
\f>\->oj 

Recalling that as \ui\ — ► 0, 

' cj° for dielectrics, 



e(r,w) 



(iw) 1 for metals 



(B6) 



[cf. Eq. l|53jl]. from Eqs. (|B4|) and (|B5|) we see that 

lim ui 2 G(r, r',u)=M, M < oo. (B7) 

|w|->0 



Needless to say, that Eq. (|B7|) is also valid for the scat- 
tering part of the Green tensor, G^(r, r', u). 



APPENDIX C: SUMMATION FORMULAS FOR 
LEGENDRE POLYNOMIALS 

The Legendre polynomials obey the relation |4jj 



J2 C nm co S (m\)P™(x)P™(y) = P n (0, 



m— 



where 



C nm — (2 S m o) 



(n— m)! 
(n+m)! ' 



£ = xy + -\/(l — £ 2 )(1 — y 2 ) cos A. 
For A = and x = y = cos 0, Eq. (|(J1|) reduces to 

n 

J2 C nmP™ (COS 6f = l. 
m=0 



(CI) 

(C2) 
(C3) 

(C4) 



Differentiating Eq. (|C1|1 twice with respect to A and 
putting A = and x = y = cos 9 afterwards yield 



z , I 



ra=0 



i(n + 1) 



(C5) 



Finally, subsequent differentiations of Eq. (|C1|I with re- 
spect to x and w and again putting A = and x = y — cos 6 
afterwards yield 



m=0 



dP™(cosl 
dfl 



n(n + 1) 



(C6) 



APPENDIX D: DERIVATION OF EQ. (fTTH) 

In Eq. (|74|l . the spherical Hankel functions h n (kor A ) 
■hn (iur A /c) can be written in the form of |40| 



i=o 



-ur A /c I 



C 

\ur A 



7+1 



(Dl) 



with some complex coefficients /ij. From inspection of 
Eqs. (Z2J and (CSJ it is seen that i?*^ and 5^ can be 
expanded in powers of u at u — 0, 



(D2) 



From Eqs. (ID1|) and (|D2|I it then follows that the inte- 
grand of the (imaginary) frequency integral in Eq. 1741) 
is a sum of terms, which are all of the same general struc- 
ture 



fjk(u) = a(iu)u 3 



uR 



c 
ur A 



fc+2 



-2uva/c 



(D3) 



(J, k are nonnegative integers). For j > (k — 1) this is 
a polynomial in u times an exponentially decaying func- 
tion. The only relevant contributions to the frequency 
integral come from the maximum of fjk(u) at a frequency 
Mo satisfying 



_d_ 

du 



fjk(u) 



0, 



thus 



u « Uq 



(j + 1 - k)c 
2r A 



(D4) 



(D5) 



where we have used the fact that a(iu) can be regarded 
as almost constant for the small frequencies considered 
here. For j < (k — 1), fjk(u) is a monotonically decreas- 
ing function, and relevant contributions to the frequency 
integral can only come from regions, where 



u < 



2r A ' 



(D6) 



because for larger frequencies the exponentially decaying 
factor becomes too small. Combining Eqs. I|D5(I and 
(|D6|I . it can be said that the relevant contributions to the 
frequency integral come from regions, where u<c/r A - In 
these regions we have 



\koR\, \kR\ ~ — < — 

C TA 



(D7) 



This means that in the long-distance limit r A R the 
main contributions to the integral come from regions, 
where |/c i?|, <C 1. We may therefore expand the 

coefficients B^' N in Eqs. J72J) and (J73J) in powers of k R, 
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on exploiting useful relations in Ref. |40| , and retain only 
the leading terms: 



B% =o[(fc 0J R) 2n+3 ] 



(D8) 



B^ ^ n ,t m ^ rT^ 1 A koR) 2n+1 - (D9) 
™ [(2n+l)!!] 2 e(iu)n + n+l K ' v ; 

In this way we find that the leading term in Eq. I|74(l 
comes from the two terms containing B^ with n = 1. 
Keeping only these terms, using 



1 



(D10) 



and changing the integration variable according to u->z 
= uta/c, we arrive at Eq. I|7t)|) . 



APPENDIX E: DERIVATION OF EQ. (l82l) 

Provided that 



n > 



4 ' 



(El) 



the spherical Bessel and Hankel functions appearing in 
Eqs. {72J ~ I7JJ can be approximated by [Zof 



and 



in (2) 



(2n+ 1)!! ' 
.(2n- 1)!! 



Z {n+X) ' 



(E2) 



(E3) 



respectively. As we know from Appendix [D] the main 
contribution to the frequency integral in Eq. I|74|) is from 
those values satisfying the condition (|D6J| . Hence, the 
condition (lElll becomes 



n > 1, 



(E4) 



because z ~ uR/c ~ uta/c in the short-distance limit. 

Substituting in Eqs. iJTJJ - J7JJ for j„(z) and /^(z) the 
expressions (|E2(1 and (|E3|) . we derive after some algebra 

B™ ~ 0, (E5) 



n(n + l)(2n + 1)-B. 



/v 



h£\k r A ) 



N 



-(2n+l)B£ 



1 



fco^A dr A 
e(m) — 1 



n 2 



n(n + l) 
(Va) 3 £H + 1 

1 e(iu)- 1 (2n + l) 2 
(k r A ) 3 e(iu) + 1 4 

1 — 1 



-2i 



n(n + 1) I — 



R 

''A 

R 

'"A 



2n+l 



2n+l 



2n+l 



.(E6) 



(/c ^a) 3 e(«u) + 1 
Whereas in the long-distance limit we could neglect all 
terms but the n — 1 one (Appendix [DJ , in the short- 
distance limit, as can be seen from Eq. (|E6|I . the param- 
eter R/r^ being very close to one, we encounter the op- 
posite extreme, where the main contribution comes from 
those terms corresponding to high orders n. The main 
contribution to the sum over n in Eq. (|74|l comes from 
the peak at m determined by 



d 

dn 



n(n + 1) ( — 



2n+l 



' n 2 e 2n(ln_R-lnr A ) 



R 



2)!i 



= 2n x [ — ) [1 + m (InR- lnr A )] = 0, (E7) 



from which we find 



ni = 



1 



R 



In r a — In i? Ata 



(E8) 



because of Ata/R <C 1. Since the main contribution 
to the sum over n in Eq. 174|) comes from values of 
n«ni3>l, where the approximate formulas (|E5|I and 
(|E6() are valid [cf. Eq. JE3J]. Therefore we introduce 
only a small error if we extrapolate these formulas to the 
terms with small n. Then the sum over n is equal to the 
second derivative with respect to (R/rA) 2 of a geomet- 
ric sum, which can be performed in closed form to yield 
Eq. (Hi- 
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